Current-induced spin polarization in nonmagnetic semiconductor junctions 
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Spontaneous spin polarization of the electrical current flowing through nonmagnetic semiconduc- 
tor junctions can be generated by carrier scattering processes that are independent of the carrier 
spin. The two required elements for current-induced spin polarization are (1) the presence of built-in 
spatially-varying electric fields in the junction and (2) energy-dependent carrier scattering processes. 
Spin-orbit interactions are not required for this effect, thus it should occur in materials like silicon 
that lack significant spin-orbit interactions. Approximate analytic expressions as well as detailed 
numerical simulations of the time-dependent nonlinear spin transport in a GaAs junction strongly 
suggest that the recent experimental observation of current-induced spin polarization in this system 
[Y. Kato, R. C. Myers, A. C. Gossard, and D. D. Awschalom, Phys. Rev. Lett. 93, 176601 (2004)] 
may be explained by this effect. 
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The goal of generating a large nonequilibrium spin po- 
larization in a nonmagnetic semiconductor has led to dra- 
matic recent advancements in understanding the funda- 
mentals of transport in magnetic and nonmagnetic semi- 
conductors, and is an essential requirement for semicon- 
ductor spintronic applications. Q, 0] Efforts to achieve 
high-efficiency spin injection from a magnetic material 
into a nonmagnetic semiconductor 0, 0, IS IE S 13 have 
yielded > 90% spin polarization from magnetic semicon- 
ductor contacts and > 70% spin polarization from fer- 
romagnetic metal contacts at room temperature. The 
recent experimental demonstration[9j, [lCJ, [llj of the spin 
Hall effect 00110 showed that spin polar- 
ization, both in density and current, can be induced in a 
nonmagnetic semiconductor through the spin-orbit inter- 
action (without recourse to any magnetic materials). The 
same spin-orbit interaction that generates these currents, 
however, also leads to rapid spin decoherence^ji so the 
spontaneous generation of spin polarization in a mate- 
rial with a very long spin coherence time remains elusive. 
One pathway towards generating spin polarization in ma- 
terials with weak spin-orbit interaction, a "spin Gunn 
effect", was suggested recently 0, but this pathway re- 
quires the semiconductor to be undergoing the rapid os- 
cillations and moving inhomogcncous domains character- 
izing electrical current flow in the Gunn effect. In addi- 
tion to the significant heating and microwave radiation 
characterizing the Gunn effect, only a very few semicon- 
ductor materials have exhibited this effect to date, and 
they all have strong spin-orbit interactions. 

Here we identify a much broader class of nonmagnetic 
materials and structures for which spontaneous spin po- 
larization will arise dynamically during carrier transport, 
but without the spontaneously-formed space-charge do- 
main of the Gunn effect. The essential requirement is the 
presence in the semiconductor of a region of static inho- 
mogeneous electric field, such as arises naturally when the 
doping of a semiconductor region is varied spatially. The 
resulting spontaneous spin polarization is unrelated to 



the spin Hall effect, as it can arise in the complete absence 
of spin-orbit interactions. Thus spontaneous spin polar- 
ization should be achievable in materials with long con- 
duction electron spin coherence times, such as is known 
for ZnO0, and is thought to be the case for silicon. 
Due to the carrier energy dependence of scattering pro- 
cesses, ubiquitous in ordinary semiconductors, a small 
excess spin polarization (thermally nucleated or seeded 
with a magnetic field) can be dramatically amplified dur- 
ing transport across the region of static electric field. 
This is demonstrated through the analytic and numeri- 
cal analysis of the nonlinear equations for time-dependent 
current transport across a junction between two regions 
of a semiconductor with different n-type doping levels (an 
n + — n junction). This effect may explain a very large, un- 
explained, experimentally-observed current-induced spin 
polarization in GaAs.|20j 

The phenomenon of spin amplification of current trans- 
port through an n + -n junction appears to be generic, and 
occurs for a wide range of temperatures and doping. Al- 
though calculations here arc presented for GaAs, similar 
calculations for other materials indicate similar spin am- 
plification (leading to larger current-induced spin polar- 
ization in materials with smaller spin-orbit interactions). 
Approximate analytic expressions for the spin amplifica- 
tion have been found which permit a qualitative under- 
standing of the temperature and doping dependence of 
the effect. Also presented here will be detailed numeri- 
cal calculations for a specific GaAs junction at 50K with 
a doping level of nci = 2.25xl0 16 cm~ 3 for x < and 
of nrj2 = 1.35xl0 16 cm -3 for x > 0. Similar numerical 
results have been obtained for many other junctions as 
well. To help the numerical stability of the calculations 
the two doping levels considered in a junction differ by 
no more than a factor of two, although the analytic re- 
sults indicate larger doping variations will be even more 
effective in generating spontaneous spin polarization. 

The larger density of carriers on the n + side of the 
junction leads to carrier diffusion across the interface and 
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Figure 1: (color online) Dependence on position of (a) carrier 
density, (b) electric field, and (c) spin polarization amplifica- 
tion rate (SPAR) for unpolarized carriers in an n + -n junction 
in GaAs at T = 50K. 



the growth of an inhomogcncous space-charge dipolc elec- 
tric field until the diffusion current is balanced by a drift 
current from the dipole field. This space-charge field will 
be the inhomogeneous electric field that drives current- 
induced spin polarization. The equilibrium charge carrier 
distribution and electrochemical potential near the inter- 
face can be determined precisely from a self-consistent 
Poisson equation, however the calculations require these 
quantities far from equilibrium where its joint solution 
with the nonlinear drift-diffusion equations would be 
required. As we must solve dynamic nonlinear drift- 
diffusion equations for spin densities, already a compu- 
tationally challenging task, we model the charge den sity 
through the junction with a simple analytic model |2lll22l | 
as shown in Fig. QJa) . The charge density is assumed to 
vary linearly on both sides of the interface. The two 
slopes and the intersection point are found by requiring 
global charge neutrality, continuity of the charge density 
across the interface, and the correct values of the electric 
fields at ±oo. This produces a parabolic analytic expres- 
sion for the electric field near the interface. Figure ^b) 
shows the calculated electric field for an asymptotic field 
of 200V/cm as x — > oo. 



The origin of current-induced spontaneous spin po- 
larization in non-magnetic semiconductors without the 
spin-orbit interaction is a dependence of the carrier mo- 
bility, jUf and /ij on the carrier spin polarization P = 
(rif — n i)/( n i + We approximate this spin polar- 

ization dependence identically to Ref. 



/i TU) = /i(n/2) • (1 + (-)aP) , 



(1) 



where the coefficient a depends on the dominant spin- 
conserving carrier scattering mechanism. A discussion 
of the dependence of a on carrier density, temperature, 
and spin-conserving scattering process can be found in 
Ref. a also depends on P, which will cause our dy- 
namical equations for spin polarization to become non- 
linear. 

To explore the consequences of the spin-polarization- 
dependent mobility in a static inhomogeneous electric 
field requires the numerical solution of spin-dependent 
drift-diffusion equations 0: 
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where T\ is the spin relaxation time of the electrons and 
D-\(l) is the diffusion constant of the electrons for spin 
up(down). Substituting /i|(x) into these from Eq. l|T]l. 
keeping only terms up to first order in P, and assuming 
D\ = Di yields 
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is the spin polarization amp lification rate (SPAR) at a 
given position. In Ref. [l8| by shifting to the moving 
frame of the Gunn domain and neglecting carrier diffu- 
sion Eq. Q was reduced to dP/dt = (7 - Tf^P, which 
was simple to solve even for a P-dependent 7. 

For an n + — n junction all terms in Eq. Q must be 
retained, which makes the calculation considerably more 
difficult, but the results much more general. For the car- 
rier density and electric field of the junction shown in 
Fig. DIa,b) the SPAR for P = is shown in Fig. QJc). 
For this junction the characteristic nature of the SPAR 
is a short region of negative values on the n + side of the 
interface followed by an extended region of positive val- 
ues on the n side of the interface. The functional form of 
the SPAR can be approximated as 



l{x) = 7max(l - x/x a ), 



(6) 
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where the length of the accumulation region (50 nm in 
Fig. ^) is x a . The polarization of a pulse of partially po- 
larized spins should be amplified when it travels through 
a region of positive SPAR (exceeding the spin relaxation 
rate). The maximum SPAR in Fig.^c) is 3 GHz, which 
leads to amplification if the spin lifetime exceeds 300 ps. 
At 50K in the GaAs-based experiments of Ref.^the spin 
lifetime is always at least an order of magnitude longer. 

For junctions between two regions of differing doping 
levels, and for the more common case0 of a > 0, each 
of n, [i, and E are decreasing in the depletion region (see 
Fig. H for the n and E dependence; for a > 0, /i follows 
the n dependence). With these trends (from Eq. ©) the 
SPAR's are positive in the depletion region. For cases 
with a < a more detailed calculation is required to 
know where the SPAR's will be positive. For ionized im- 
purity and optical phonon scattering, the common mech- 
anisms of mobility in semiconductors at low and high 
temperature, a > 0. The position dependence of the 
SPAR is similar for all such junctions, an initial suppres- 
sion region followed by an extended region of positive 
values, linearly diminishing up to the edge of the accu- 
mulation region. 

The behavior of a wide variety of junctions can be es- 
timated from the maximum value of the SPAR, which is 
shown in Fig. El AH the n + — n junctions considered have 
?i£>2 = 0.6ri£)i and the mobility is considered to originate 
from ionized impurity scattering. As the temperature is 
increased the dependence of the mobility on spin dimin- 
ishes (a decreases ^^) so the SPAR diminishes as well 
(Fig-GI a ))- As the carrier density is increased the depen- 
dence of the mobility on spin increases (a increases [lg) 
so the SPAR increases (Fig.^b)). At room temperature 
in a system doped n = 10 18 cm~ 3 the SPAR exceeds the 
spin relaxation rate if T\ > 5ps. 

For the polarization of a pulse of carriers to be signif- 
icantly amplified in the junction the pulse must remain 
in the region of positive SPAR for sufficient time. The 
instantaneous velocity of a piece of a pulse of carriers, 
from the coefficient of the dPj dx term of Eq. (@J is 
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For the junction considered in Fig. ^ti(j;) < fiE in the 
region of amplification, < x < 50 nm, both because 
a > and dn/ dx < 0. This self-consistent reduction in 
the effective velocity of a carrier pulse plays an impor- 
tant role in amplification, as otherwise the pulse would 
only remain ~ 25 ps in the amplification region. An an- 
alytic expression can be derived with the assumptions of 
Fig. Hand (1) the SPAR has the form of Eq. ©, (2) the 
polarization dependence of a (and thus of 7) can be ne- 
glected, and (3) diffusion (the d 2 P/dx 2 term in Eq. (@J 
can be neglected. Under these conditions v(x) = i> e ff is a 
constant in the positive SPAR region. The polarization 
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Figure 2: (color online) The maximum spin-polarization am- 
plification rate (SPAR) 7 ma x, for unpolarized carriers as a 
function of temperature and junction density. Ionized impu- 
rity scattering dominates the mobility. 



of the pulse would be 

P{t) = Pinitiaiexp [( 7max - T^)(t - v e gt 2 /2x a )] . (8) 
Evaluating at t — x a /v c ff yields an amplification of 
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A numerical calculation that includes the spin polariza- 
tion dependence of the SPAR and v e g, and also includes 
diffusion, for the junction in Fig. is shown in Fig. |3| 
at T = 50K. The system is initialized at t = with a 
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weakly spin polarized gaussian pulse centered at x = —d\ 
on the n + side of the n + —n junction. Subsequent panels 
in Fig. are labeled by the time and indicate that the 
pulse is suppressed as it propagates through the deple- 
tion region (x < 0) , as expected from the negative SPAR 
in Fig. ^c) in this region. During this the pulse drifts to- 
wards the interface and spreads due to diffusion and the 
position-dependent velocity v(x). At t = 32ps the spin 
polarization reaches the amplification region at x > 
and begins to increase. By t = 160ps the spin polariza- 
tion has saturated at P ~ 1. The distorted pulse then 
slowly drifts towards the edge of the accumulation region. 
Outside of the accumulation region the velocity is much 
larger than within, so it takes ~ 5ns for the highly spin 
polarized pulse within the accumulation region to drive 
a highly spin-polarized pulse out of the junction. Then 
this spin polarization pulse packet moves away from the 
junction. The output appears to be a nearly 100% spin- 
polarized pulse that is several ns wide (corresponding to 
~ 10/im in spatial width). 
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Figure 3: Propagation of an initial weakly spin-polarized 
pulse through the n + — n junction of Fig. 1. Amplification 
occurs for x > and yields a nearly 100% spin-polarized pulse 
into the n-doped region. Fhe spatial dependence of the spin 
polarization is shown for 10 time slices from 8 ps to 16 ns. 

We found the spontaneous generation of highly spin- 
polarized pulses of carriers occurs in current flow through 
an n + — n junction. An analytic expression was derived 
for such junctions, neglecting diffusion and spin polar- 
ization saturation. Numerical results also support these 



findings. Calculations showing the presence of highly 
spin-polarized current in GaAs junctions indicate these 
results provide a possible explanation for the current- 
induced spin polarization of Ref. |U 
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